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Using NLL values for Wilson coefficients and including the contributions from the penguin 
diagrams, we estimate the amount of nonfactorization in two-body hadronic B decays. Also, we 
investigate the model dependence of the nonfactorization parameters by performing the calculation 
using different models for the form factors. The results support the universality of nonfactorizable 
contributions in both Cabibbo-favored and Cabibbo-suppressed B decays. 
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I. INTRODUCTION 



S3 . 

The decay rates (and, in few cases, polarization) of a large number of hadronic channels have been experimentally 
measured to sufficient accuracy. However, as yet there is no reliable theoretical method to derive the corresponding 
amplitudes starting from the basic principles of the Standard Model. This is due to our inability to quantify the 
participation of strong interactions in such processes. The powerful theoretical tools of perturbation theory which are 
used in purely electroweak interactions are not very useful in situations involving strong interactions. Nonperturbative 
techniques, such as lattice calculations and QCD sum rules, are still under development. However, the asymptotic 
freedom property of QCD allows us to separate the gluon contribution, in a given process, into that due to high 
energy (hard) gluons and that due to low energy (soft) gluons. The former contribution is relatively easy to compute 
using perturbation techniques and renormalization group equations. In fact, an impressive amount of work in this 
regard has been done M-n , where hard gluon effects were parametrized through Wilson coefficients which have been 
0\ ■ calculated up to next-to- leading logarithmic order (NLL) . It is the soft gluon contribution which is difficult to handle 
0\ | and constitutes the main source of uncertainty in hadronic weak decays. 

In general, the effective Hamiltonian in the hadronic decays of B mesons takes the form H e g ~ ^ CiQi where Cs 

O ,' are the Wilson coefficients that contain the hard gluon (or short distance) effects and Qs are four-quark operators 

JL • that are products of two Dirac currents. In calculating the decay amplitudes for two-body hadronic decays we then 

encounter matrix elements of the form (/i f 2 \Q\'i), where i is the initial state particle, f\ and f 2 are the final-state 

-— ■ particles. For example, the effective Hamiltonian relevant to the process B — > D + ir~ is given by (see the next section 



for details), 



H ctf = %V cb v: d [d Qi + C 2 Q 2 ] , (1) 



where 

Qi = (ci bi) L (dj Uj) L , (2) 

Q 2 = (cib^L (djUi) L 

= jr(ci bi)L (dj u 3 )l + ^(ci^i]bj)L (d~ k \ kl ui) L . (3) 

The subscripts i, j, fc, and I represent color, N c is the number of colors, As are the Gell-Mann matrices and L represents 
a left-handed current, 

(016^ = 5^(1-75)^. (4) 

In a diagrammatic language similar to that used in Ref . M , we show in Fig. p] different diagrams contributing to the 
decay amplitude of this process. Regarding these diagrams several remarks can be made. First, the leading diagram 
generated by Q\ (Fig. [j] (a)) is factorizable, where all soft gluon exchanges inside the closed loops are accounted 
for by the B — ► D form factors and the pion decay constant. Second, the next-to-leading diagram generated by ^ 
(Fig. [l] (b)) is not factorizable and is of the order of l/N% relative to the leading one. This is due to the double soft- 
gluon exchange between two color-singlet mesons. Third, when the operator Q 2 is rewritten using Fierz transformation 



and color algebra in the form of Eq. (||) , it generates, the lower set of diagrams in Fig. ffl. The diagram in Fig. [I] (c) 
which is generated by the color-singlet current is factorizable whereas those in Figs, [l] (d) andjl] (e) are not. Finally, 
the annihilation diagram in Fig. g is generated by Q 2 alone. 

As we said earlier the factorizable part can be parametrized in terms of form factors and decay constants. On the 
other hand, the nonfactorizable parts are more difficult to handle. A popular method is to omit the nonfactorizablc 
parts altogether and to assume complete factorization of the decay amplitude. Some justification for this approach 
is provided by the fact that the nonfactorizable diagrams are nonleading in the 1/N C expansion. However, there are 
two weak points in this argument. First, 1/N C (where N c = 3) is not a good expansion parameter. Second, even the 
factorizable part carries some uncertainties due to the model dependence of the form factors and scale dependence of 
the Wilson coefficients. Of course, experiment has the final say regarding the validity of such an approach. 

In B decays into two psedoscalar mesons, the annihilation diagram contribution can be argued to be small. For 

example, in the process B — » D + ir~ we expect the ratio of the annihilation diagram to the factorizable tree diagram 
to be of the order (here the annihilation diagram is also calculated in the factorization approximation) 

C 2 (ml -ml) f B F D ^(m B ) _ f B F D ^(m%) 

~ d (m% -ml) U F B ^K) ~ ' U F B ^ D (m%Y K } 

where the values of the Wilson coefficients used are shown in Table |. So, unless F "* 71 (m B ) has a very high value 
(which is not expected) the contribution of the annihilation diagram is, in all likelihood, very small. In this paper, 
we neglect the annihilation terms in decays involving vector particles also. 

At this point we should differentiate between two approaches to factorization. In the first approach, called naive 
factorization (see for example Ref. [g|), both leading and sublcading diagrams are included in amplitude calculations. 
In the second approach (see for example Ref. 0), called the large iV c limit, only the leading diagrams in l/N c are 
included. If the leading diagram in a decay process is generated by the operator Q\, factorization works at its best. 
This is because the subleading diagrams are suppressed further by the Wilson coefficient Ci. In the large N c limit, 
the decay amplitude for such processes, which are labelled class I M, is proportional to 

CxWKXMJ^O) (6) 

In this equation J^ and J' are two Dirac currents. In naive factorization the amplitude is proportional to 

Ci + ^r) WKX/al-W- ( 7 ) 

The difference between decay rates predicted by these two amplitudes (defined by Eqs. (||) and (M))is not very large. 
It is about 20% (15%) in the case of D (B) decays. However, in class II processes, where Q2 generates the leading 
factorizable diagram, we find that the decay amplitudes are proportional to 

C 2 (f2\J^)(fi\J',\0) (8) 

and 

(d + ^) WliX/iW (9) 

in the large N c limit and in the naive factorization, respectively. Unlike processes of class I, the inclusion of the 
subleading factorizable diagram in class II processes has a significant effect on the decay amplitude. For example, in 
the case of D decays the predicted decay rates drop by one order of magnitude. In class III processes, both Q\ and 
Q2 generate leading and subleading factorizable diagrams. So, the amplitude is proportional to 



and 



c 1 </ 1 |j^N)(/ 2 |j;|o) + c 2 </ 2 |^| 1 )(/ 1 |j;|o) (10) 

(ci + ^) (/iTOWlo) + (c 2 + ^ (/ 2 |J"|»X/il4|o) (if) 

in the large N c limit and in the naive factorization, respectively. 

When factorization was applied in studying the decays of D mesons, it was found that the large N c limit worked 
much better than naive factorization. This seemed logical from the point of view of 1/N C expansion because J7|] it is not 



justified to include only part of the terms in a certain order. This led to the anticipation that 1/N C expansion can also 
be used to describe B decays. However, when B measurements started to appear it was found that both factorization 
approaches failed to describe the new experimental results, especially for class II and class III processes. For example, 
the decay rate of the process B — > KJ/ijj was found to be about one order of magnitude higher than that predicted 
by factorization |0,M. This shifted the attention of many physicists towards the nonfactorizable terms, which were 
previously omitted. The argument used is that these terms are not negligible and they contribute significantly to the 
decay amplitude. However, in the case of D decays it just happened that the subleading factorizable diagram almost 
canceled out the nonfactorizable one due to a destructive interference between the two. On the other hand, in the 
case of B decays the interference seemed to be constructive. 

Ideally, if the nonfactorizable contributions could be calculated from the basic principles of the Standard Model we 
could systematically test the validity of the 1/N C expansion. Unfortunately, with the present theoretical capabilities 
this is a difficult task (see Refs. [[tl],[l2| for some attempts in this direction). However, until such computational methods 
become more reliable, we can still do phenomenological parameterizations that can be compared with experimental 
data. 

So, in going one step beyond factorization we can write the decay amplitude up to the first order in 1/N C but includ- 
ing the nonfactorizable diagram generated by color-octet currents. Changes would then be made to the amplitudes 
m &, (H) an d C P-iP through the replacements: 

C a + %f) -^(c 2 + ^ + C 1 e s ). (12) 



NJ V N, 

where £% represents the contribution of the color octet nonfactorizable diagram relative to the factorizable one. In 
Refs. [|l3H16| several authors attempted to estimate the amount of nonfactorizable contributions in several decay 
channels of D and B decays. The results highlight two main points. First, in D decays the value of £s is negative 
and highly process dependent. Second, in B decays eg is positive and is not as process dependent as in D decays. In 
addition, the results obtained in Ref. pl| seems to satisfy the inequality 

\e s (D -^VP)\> \s 8 (D -^PP)\> \e s (B decays) | , (13) 

which finds some justification in the color transparency argument Jl7| ], According to this argument soft-gluon exchange 
becomes less important as the final products of an interaction become more energetic. This is because the final state 
particles leave the interaction region very quickly allowing little time for final state interactions. It should be mentioned 
here that some authors |18|-|20| treated N c as a free parameter whose effective value is used to indicate the amount of 
nonfactorization in a decay process. 

One of the issues that arises when parameterizing nonfactorizable contributions is how to handle processes with two 
vector mesons in the final state? This is because it is not clear whether or not the three Lorentz scalar structures of 
the decay amplitude should receive the same contribution from the nonfactorizable terms, i.e. does nonfactorization 
lead verily to an over all factor? This issue was discussed in Refs. 14 13[|. In Ref. |2l| we tackled this issue in some 



detail for the process B — v J/ipK*. Using a full amplitude measurement |22| by CLEO, the amount of nonfactorizable 
contribution to each of the three Lorentz-scalar structures was calculated in five different models for the form factors. 
The results allowed an explanation of the experimental data using equal amount of nonfactorization in each part of 
the Lorentz amplitude, implying that an overall nonfactorization factor was adequate. 

Assuming universality (process-independence) of the nonfactorization parameters in B decays, we estimated their 
values in Ref. [£3[ using a more definitive calculation. However, for the Wilson coefficients we used the values 
calculated up to leading logarithmic order (LL) and neglected all contributions from the penguin diagrams. Also, 
all the calculations were done using only one model for the form factors (BSW II model). The results supported 
the proposition of universality of the nonfactorization parameters in Cabibbo-favored B decays. However, a number 
of questions were also raised. First, how will the results change if we use the Wilson coefficients calculated up to 
NLL? Second, how do we parameterize nonfactorization generated by penguin diagrams and how important they are? 
Third, how much model-dependence is there in the estimated nonfactorization parameters? Finally, can we extend 
the proposed universality of the nonfactorization parameters to include Cabibbo-suppressed B decays? These are the 
questions we try to address in this paper. 

The paper is arranged as follows in Sections II and III we present the Wilson coefficients and CKM matrix elements 
used in the calculations. In Section IV we calculate the effects of penguin diagrams and NLL Wilson coefficients 
on the predictions of naive factorization. In Section V the nonfactorization parameters are estimated in five models 
for the form factors. In Section VI we show branching ratio predictions of several sets of Cabibbo-favored and 
Cabibbo-suppressed B decays. The last Section is a discussion of the results and a conclusion. 



II. WILSON COEFFICIENTS IN NLL 

In the absence of strong interactions, the effective Hamiltonian for the process b — > ccs, is given by 

G F 



(14) 



When QCD effects are included, the contribution of the penguin diagram in Fig. p| should be considered beside the 
current x current diagrams. As a result, the effective Hamiltonian generalizes to |lq,E(| 



W e ff 



Gf 

7f 



v ab v: 3 (Ci q\ + c 2 Ql) + v cb v; s (d Q\ + c 2 qd + (v ub v; s + v cb v; s ) J2 c % q, 



where 



Qi = {uibi) L (sjUj)l 
Q l 2 = {uibj) L (sjUi) L 
Qi = (cibi) L (sjCj)l 
Q c 2 = (cibj) L {sjCi)L 
Q-3 = (s l b l ) L ^2(qjqj)L 

Q 

Qa = (s. t bj) L ^2(q J q l )L 
1 

Q.5 = (sibi) L ^2(mj)R 

q 

Qe = {Sibj) L 22(qjQi)R- 



i=3 



(15) 



(16) 



The superscript R represents a right-handed current such that 

fe«i)fl =«j7 /1 (l + 75)«j- 



(17) 



Even though, the local operators Q\ and Q 2 don't contribute to processes of the type b — ► ccs through tree diagrams, 
they do contribute through penguin diagrams. In writing (flq), we have made use of the following unitarity condition 
of the CKM matrix elements 



v ub v* a + v cb v; s + v tb v t * s = o. 



(18) 



The Wilson Coefficients, C\ and C 2 , are determined by requiring the effective Hamiltonian to reproduce the am- 
plitude calculated in the full theory. In LL calculations this matching of the full theory onto the effective theory is 
done by keeping only the terms containing logarithmic corrections of order a s ln(). In NLL calculations, the terms 
containing constant contributions of order a s are kept in addition to the terms of order a s ln(). The Wilson coefficients 
turn out to be regularization scheme dependent. However, if the matrix elements (Qi(fJ.)} are evaluated at the same 
scheme as the Wilson coefficients Ci(fi), the scheme dependence cancels out. From Eq. (|15|), we see that the decay 
amplitude in the effective theory has the form 



As oc Ci(fi) (Qi{fj,)) ■ 
If these matrix elements are related to the tree matrix elements by 

(Q»(m)) =9ij(v) (Qj) tICC ■ 

then we can write the effective amplitude as 

Acs oc Ci(p) gij(p) (Qj) t! 
oc Cf (Q^ 00 . 



(19) 
(20) 

(21) 



At the quark level, the scale and scheme dependence of (Qi(n)), which is carried by g(/J,), cancel the scale and scheme 
dependence of the Wilson coefficients [[18 20 2^]|(|. So, both Cf s and (Qi) ree are scale and scheme independent. 
From above, we can write the decay amplitude as 



(H cS ) -J 



v cb v; s (cf{Qir ee + cf(QT' cc ) + E (*W. cr cff + v cb v; s cr a ) (Qi 



i = 3 



(22) 



where the penguin contributions from (Q" =1 2 ) and (Q| =1 2 ) are included in Cf^ 6 and Qiff 6 , respectively. 
Using the conventions of Ref . |18J , the explicit forms of C° ff are given by 



Cf 1 = Ci 



C| ff = C* 2 



3 — G 3 



z^g off ^ 

o 4 — O4 



a 



q cff 



c 3 



4tt 
24tt 



C, 



ry + 7y In — , „ 3 . 
f* J ij 

r v + ~{ v In I 



C, 



/" / 



2j 



[C t (m g ) + C p ] + 



47T 



7 v ln 



M / 



8?r 
24tt 



C 7 



•■y 



P.K) + c p ] + ^( rv : + ri in ^) c 

47r V M /4j 



[C t (m g ) + C* p ] 



47T 



7 v ln 



m b \ 

m / 



c. 



Sj 



c; 



q cff 



C4 + ^[C t (m g ) + C p ] 

87T 



4tt 



7 y ln 



r;? 5 



" ' 6j 



c,- , 



(23) 



where the superscript (7 represent either the u or the c quark. In (|23|), the quantities Ct(m q ) and C p arise from the 
penguin-like diagrams of the operators Qf =1 2 an d Qi=3,...,6, respectively, whereas the matrix ( ry + 7y In ^ J is due 
to the vertex and self-energy corrections to the operators Q»=i,. e- These quantities are given by [Q 



(I 



rv 



-7 \ 

-7 I 

I -7 

-7 I 




V 



— I 1 

-3 f 



(24) 



lv 



(- 2 
6 







Vo 



6 

-2 











-2 
G 
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-6 



(25) 



-16/ 



C t {m q )=C 1 (ji) 



C P = C 3 (n) 



2 2 , m^ A / fc 2 

3 + 3 ln #" AFl (^ 



2 to 
In 



3 /i 
+ [C 4 (m) + C 6 ( A )] 



2 TO 
- m — ; 



2=U,d,S,C,fr L 



2 



2 TO, 

-In— |- 



AJi -* 



AFif 


4)1 


V 


'fc 2 y 




,™f/_ 


1 



(26) 



(27) 



and 



-4 / dx x(l — x) ln[l — z x(l 
Jo 



AFiO) = -4/ dx x(l-x)\n[l- zx(l-x)}. (28) 

Jo 

Here, k is the momentum carried by the gluon in Fig. H. 

The calculation of the penguin-driven amplitudes in the factorization assumption involves additional assumptions, 
an effective value of k 2 for example. In a complete calculation |27| k 2 would not be a variable; it would be integrated 
over the wave functions of the hadrons with its own uncertainties. In the absence of a complete knowledge of the 
hadronic wave functions, the choice is either to select k 2 judiciously or to admit new unknowns through the hadronic 
wave functions. In penguin calculations, one generally opts for the first alternative and chooses k 2 in the range: 
TO &/4 < k 2 < m 2 /2. In the calculations presented here we have chosen k 2 = rn 2 /2. 

Working in the naive dimensional regularization (NDR) scheme, the Wilson coefficients at the scale p = 4.6 GeV 
are listed in Table |. Using the definitions in Eq. (|2J|) and the values of Ci (/i) in Table | we calculate in Table [0] the 
values of the effective Wilson coefficients C° ff . For C^ s — Cg ff , which include the contributions from the QCD penguin 
diagrams, we list two sets of values, depending on which quark is in the loop (see Fig. [|). For quark masses, we used 
the the following running values at the 6-quark mass scale [p0[: ra u — 3.17 MeV, nid = 6.37 MeV, m s = 0.127 GeV, 
m c = 0.949 GeV, m b = 4.34 GeV and m t = 170 GeV. 

III. CKM MATRIX ELEMENTS 

To be able to calculate the decay amplitudes of the processes considered in this work, we need the values of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. In Wolfenstein parametrization pqj, and up to the fourth 
order in A, these elements take the form 

V ud V us V ub \ ( 1-±A 2 A A\ 3 (p-ir)) 

V cd V cs V cb \ = \ -X{l + iA 2 X i r 1 ) 1-|A 2 AX 2 | . (29) 

V td V u V tb I \A\ 3 (l-p-ir}) -AX 2 1 



In a recent update |29(| , the four parameters in Eq. (|29| ) were given the values A = 0.2205, A = 0.81, p = 0.05, and 
7] = 0.36 using a x 2 fit to experimental data. The needed CKM matrix elements would then be: 

V ud = 0.976, V us = 0.221, V ub = 0.00316 e~ 1A3 \ 

V cd = -0.221, V cs = 0.976, V cb = 0.0394. (30) 



IV. DECAY RATES IN NAIVE FACTORIZATION 
A. Type b —* cud Decays 

Decays of type b — > cud proceed via the effective Hamiltonian presented in Eq. (|l|), where only tree diagrams 

contribute to the decay amplitudes. In processes of class I (for example B — > D + ir~) the factorizablc part of the 
amplitude is proportional to ai = C\ + C2/3. As can bee seen from Tables | and ||, the value of this parameter in 
LL (a\ h = 1.032) is almost the same as that in NLL (a^ LL = 1.036). Since the decay rates are proportional to |oi| , 
the above values for a\ gives a difference (Ai3f ac ) of less than 1% between the branching ratios calculated using the 
Wilson coefficients in LL and the branching ratios calculated using the Wilson coefficients in NLL (see Table |lj). 
In processes of class II (for example B — > D°ir°), the factorizable part of the decay amplitude is proportional to 
<22 = Ci/3 + Ci- This parameter takes the values ( a? L = 0.090) and (a^ LL = 0.059) in LL and NLL, respectively. As 



i2_ 
a result, the predicted branching ratios (see Table III) in the naive factorization approximation, drop by about 57% 



when working in NLL. 

The decay amplitudes for class III processes (for example B~ —> D°ir~) receive contributions from two tree diagrams 
causing a dependence on both a\ and ai. Therefore, A2?f ac varies from one process to another in this class. However, 



as can be seen from Table III these changes are relatively small (less that 6%). This is caused by the dominance of 



the part of the amplitude proportional to a\ over that proportional to 02. It should be mentioned here that, unlike 



the other two classes, A£>f ac in class III processes is model dependent and the values presented in Table III were 
calculated based on the BSW II model, to be introduced later in this paper. However, due to the dominance of one 
part of the decay amplitude over the other, this model dependence is not very strong. 



B. Type b — ► ccs Decays 



Decays of type b — ► ccs, receive contributions from both tree and penguin diagrams and the relevant effective 
Hamiltonian is given by Eq. (|l5|). Consider the following four class I decay channels, B — > D + D~ , D + D* s ~ , D* + Dj 
and D* + D*~ . In order to extract the factorizable contributions to these processes, an appropriate transformation is 
needed for the operators Q§ and Qr containing right-handed currents. Starting by Q5, we can use Fierz transformation 
and color algebra to rewrite it as JL8| 



Qb = (sb) L (cc)r 

= [SiTfitO- ~ y&)bi] [Ci7"(l + Js)cj] 
= -2^(1-75)64] [5,(1 + 75)^-] 

- A [2(1 - 75)6] [«(1 + 7s)c] - [cA Q (l - 75 )&] [sA a (l + 7s )c] . 



(31) 



For Q61 we just use Fierz transformation to write it in the form 



h = (sibj)L {c 3 Ci) R 

= [s^t 1 - 75)&j] [c^C 1 + 7b)ci] 
= -2 [c(l - 75)6] [s(l + 75 )c]. 



(32) 



However, from Dirac Equation we can easily derive the following relations between (S + P) and (S — P) Dirac bilinears, 
and V and A bilinears 



[s(l + 75)c] = i 
[c(l - 75)6] = i 



9 M (s7 M c) 9 M (s7 M 7 5 c) 



m c — m s m c + m s 



m b — m c m b + ti c 



(33) 

(34) 



Using the relations ( pl[]34| ), we can write (assuming naive factorization, symbolized by the subscript fac) the decay 
amplitude for the processes mentioned above as 



A iac (B°^D + D;) = ^ 



E 



m 



D, 



V cb V; s a 1+ ^ V qb V* s [a\ + 2al 

7 -^ ' [mb — m c )(rn c + m s ) 



xp+K^iF-XDri^cKio), 



(35) 



Aac(B°-^ + ^D = ^| 



V<*t£ ai + E *W- a 4 + 2a > 



o D , 



q=u,c 



{m b - m c )(m c - m s ) 



i-W) 



x (U+K^ilB") (D* s ~ \(sc) L \0), 



(36) 



Aac(S°-^ + ^) = ^ 



W>i + E V qb V q * s (al-2a, 



g— u,c 



(m b +m c )(m c + m s ) 



x(D* + \(cb) L \B)(D-\(sc) L \0) 



(37) 



and 



A iac (B a -> d* + d; 



Gf_ 

V2 



V cb V c * s a 1+ Y, V qb V;Aat + 2a. 



• d : 



q=u,c 



(m b - m c )(m c - m s ) 



x(D*+\c 1 ,b\B)(D;\(sc) L \Q) 



Gf 

T2 



V cb V c * sai + ]T V qb V* s (al-2a 



q=u,c 



(m b +m c )(m c - m s ) 



x (D*+\ry^ 5 b\B )(D;\(sc) 1 



(38) 



where 



4 = ci 



,« _ 



,1 - 



ci 



ci 



a\ — (7 6 



o 4 , 



\C 



3) 






u,c. 



(39) 



Note that N c — 3 is used for the penguin amplitudes also. 

If the penguin-generated terms were omitted, the effect of working in NLL instead of LL would be very small. In 
fact, ASf ac is less than 1%, the same as that calculated above for the color-favored decays of type b — ► cud. However, 
if the contributions from the penguin diagrams are considered we get relatively large effects. The branching ratios for 
the processes B — ► D + D~ and B — > D + D* s ~ get reduced by 27% and 36%, respectively, while the branching ratios 
for B — > D*D S and B — » D*D* are increased by 6% and 15%, respectively (see Table |l|). To demonstrate the cause 
of these large changes, let us consider the decay B — > D + D~ . In a rough calculation, we substitute the following 



approximations in (35): a\ rs 1 



m c )(m c 



1. The change 



-0.04, ag « -0.05, KbK** « and m^ s /(m h 
in the amplitude due to penguins is then about (1 — 0.04 — 2 x 0.05) 2 -1» —27%. 

The decay amplitudes of class II processes are simpler to derive. For example, the factorized decay amplitude for 
the process B — ► K J/t/j is given by 



•Af ac (B 



K J/ijj) 



V2 



V cb V* a 2 



E v * v - 



qa 



(«! + «!) 



(K°\(sb) L \B°)(J/i;\(cc) L \0). 



(40) 



Other processes of this class have similar amplitude forms. When the Wilson coefficients generated by the penguin 
diagrams are set to zero, the branching ratios calculated in NLL are 57% lower than those calculated in LL (similar 
to class II processes of type b — > cud). The penguin effects, however, turn out to be very small (about 0.2%). This is 
because the values for 03 and as are very close in magnitude and have opposite signs resulting in a mutual cancellation. 

C. Type b — » ccd Decays 

The effective Hamiltonian for processes of type b — > ccd is similar to Eq. pq ) except that the s flavor is replaced by 
the d flavor. So, the effective Wilson coefficients are calculated from Eqs. (123), using the following expression for C p 



C P = C 3 (fi) 
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i—u.d,s,c,b 
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2 -m 
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ml 


[1- 
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1^~ 


AFi 


(- 



AFi 



AFi 



[C 4 (/z) + C 6 ( M )] 



(41) 



The numerical values for C^fj 6 turn out to be the same as those in Table O without noticeable changes. The decay 
amplitudes for the processes B — > DD~ , DD*~ , D*D~ and D*D*~ can be written from Eqs. (p5|)-(p8|) by replacing 
the s flavor by the d flavor. By doing a similar replacement, the amplitude for the process B — ► TtJ/ip can be written 
from Eq. (EoJ) . In Table III , we show the NLL and Penguin effects on the calculated branching ratios of these two 
sets of processes. 



V. ESTIMATION OF NONFACTORIZATION PARAMETERS 



In Ref. |16| , nonfactorization was parametrized through e\ and eg. These two parameters represent the size of 
the color-singlet and color-octet nonfactorizable diagrams relative to the factorizable one (see Fig. ph. By assuming 
universality (process-independence) of these two parameters, we estimated their values in Ref. P3| for Cabibbo-favored 
B decays. The estimate was done using the available experimental branching ratios for two sets of class I processes 
(B° -> D+TT-,D + p-,D + a^,D*+TT-,D* + p-,D* + a^ and B -> DD 8 ,DD*,D*D 8 ,D*D*) and one set of class II 



processes (B — ► KJ/ip, Kip(2S), J/ip, K*ip(2S)). For the Wilson coefficients we [£3| used the values calculated up to 
LL and neglected all contributions from the penguin diagrams. Regarding the form factors, we used the predictions of 
BSW II model. In general, the estimated values of E\ and e% improved the agreement with experimental measurements 
when the notion of factorization was extended to include other channels of B and B s decays. This supported the 
assumption of universality of these parameters in Cabibbo-favored B decays. 

According to Eq. (12), nonfactorization contribute to the decay amplitudes of class I and class II processes through 
the multiplicative facors (omitting the penguin contributions for simplicity of argument) 



ii 



l + £l 



9i 



£8 



l + £l 



1 



and 



C 

6 = ( 1 + e a + — e 8 



£8 



20 e 8 



(42) 



(43) 



respectively. The difference between these two factors is in the coefficient of the color-octet parameter (see Table ||) . 
Consequently, in case of class II processes the long distance effects in eg are greatly enhanced by the short distance 
effects arising from C\/a2- This is in addition to the enhancement of e% over £i, by a factor of N c , according to the 
rules of QCD. As a result, it may be harmless to ignore the contribution due to the color-singlet parameter e\. On the 
other hand, in class I processes we notice that £s is suppressed by approximately a factor of 1/3 due to short-distance 
effects. Since this compensates for the enhancement due to QCD, it is not justified to omit e± in this class and the 
two parameters should be treated on equal footing. From the preceding discussion we infer that class II processes are 
sensitive probes of Eg, whereas £\ is mainly determined by class I processes, albeit with less sensitivity. 

In this work, we re-estimate the values of the nonfactorization parameters E\ and £s using a \ 2 fit to the experimental 
branching ratios of the three sets of processes mentioned above. However, for the Wilson coefficients we use the values 
calculated up to NLL and include, in the color-singlet part of the amplitude, the contributions from the penguin 
diagrams. As for the color-octet part, we include the contributions from the operators Q\ and Qi only. For example, 
the decay amplitude for the process B — > KJ/ip, after including the nonfactorizable contributions will read 



A(B -» KJ/ip) 



Gf_ 

V2 



VcV* 



«2 



E 



VkV* [a3 

v cbv cs 




C\ £8 



(K\(sb) L \B)(J/iP\(cc) L \Q). (44) 



Since the nonfactorization parameters are model dependent, we repeat the % 2 fit using five different models for the 
form factors. The first of these models is the original Bauer, Stech and Wirbel model (called BSW I here) where 
the form factors are calculated at zero momentum transfer and extrapolated using a monopole form for all the form 
factors. The second model (called BSW II here) differs from the first one by using a dipole form to extrapolate the 
form factors F±, Aq, Ai and V. This is motivated by the consistency relations in the infinite quark mass limit derived 
in Ref. [p2[ . In the third model 32 (called NRSX here) we use the heavy-quark effective theory predictions for the 
heavy-to-heavy form factors. For the heavy-to-light form factors we use the same values as those predicted by BSW II 
model. In the fourth model (AW), developed by Altomari and Wolfenstein [[33|, the form factors are evaluated at 
the zero-recoil point corresponding to the maximum momentum transfer and then extrapolated down to the required 
momentum using a monopole form. In the last model (ISGW), by Isgur, Scora, Grinstein and Wise p4| , the form 
factors are calculated at the maximum momentum transfer and extrapolated down with an exponential form. See 
Table IV for a comparison between the form factors predicted by these models for the process B — > DD S . 

Working, for example, in BSW II model we show in Fig. (a) a contour plot of % 2 in £i-£s space. The four minima, 
labeled 1, 2, 3 and 4, that appear in this figure is to be compared with the corresponding regions in Fig. 2 of Ref. [23|. 
The value of % 2 per degree of freedom (% 2 /d) for these minima is 0.6. In an argument similar to that used in Ref. [23|, 
we exclude solutions 1, and 2 due to the severe violation of the approximate relation E\/e% = ±1/N C suggested by 
1/N C expansion. Also, solution 4 is excluded because it produces a negative value for the effective 02 parameter. So, 
we end up with the estimate E\ — —0.040 ± 0.024 and £8 = 0.137 ± 0.006. The uncertainties correspond to A% 2 = 1. 
The predictions of the other models are shown in Table [V|. 

Out of the five models considered we note that three (BSW I, ISGW and AW) do not produce a good fit to the 
experimental data. This is indicated by the high values of \ 2 /d. The remaining two models, on the other hand, 
produce a good fit to the data with x 2 /d equal to 0.6 for BSW II model and 0.4 for NRSX model. In these two models 
which show interesting fits the heavy-to-light form factors are calculated in the same way. So, the set of class II, 
processes (B — ► KJ/ip . . . etc.) have the same values for the form factors. This is reflected in close predictions 
of £3 by both models, which is not the case for the other parameter where the two sets of class I processes take 



different values for the form factors in the two models. Another point to be noticed is that NRSX model predicts 
a destructive interference between the color-singlet and color-octet nonfactorization contributions causing almost a 
complete cancellation between the two for class I processes. This is not the case for the other models which suggest 
a constructive interference. 

VI. PREDICTED DECAY RATES INCLUDING NONFACTORIZABLE CONTRIBUTIONS 

Using the values estimated for the nonfactorization parameters in a scheme with NLL Wilson coefficients and 
penguin contributions, we calculate the branching ratios, in BSW II model, for all processes considered in the previous 



work p3J. The results are shown in Tables VI and VII. Also, by assuming that the universality of E\ and eg extends 



to Cabibbo-suppressed processes, we evaluate the branching ratios for a set of class I proc esses of type b — > ccd and 



for another set of class II processes of the same type. The results are shown in Tabic |VII|. These two sets were not 
considered in the previous work. 

Our predictions show a good agreement with available experimental data which includes the recently measured 



decay channels ~W —* D* + D*~ 35 and B~ — > Tr~J/?p [|3l). In the calculations, the two states \rj) and \rf) are treated 
in the same way as in plfl where the mixing angle and wavefunction normalizations are properly taken care of. As 
for the decay constants we used the values adopted in Ref . [£3| . 

VII. DISCUSSION AND CONCLUSION 

In Ref. |23|| , we demonstrated that naive factorization, in LL and using tree diagrams only, gives reasonable pre- 
dictions in comparison with experimental measurements for the branching ratios of class I and class III processes. 
However, for class II processes the predicted branching ratios were very low. Also, it was demonstrated that including 
nonfactorizable contributions through the parameters £\ and eg improves considerably the predicted branching ratios 
for the latter class while preserving the reasonable predictions for the other two. 

In this work, we find that by working in NLL (with no penguins) the predicted branching ratios (in naive factor- 
ization) of class I and class III processes are very close to the LL predictions and to the experimental values (see 



Table III). For class II processes, on the other hand, the predicted branching ratios are considerably lower (by about 
57%) than the LL predictions making the disagreement with the experimental values even worse. However, this prob- 
lem is greatly remedied by including nonfactorizable contributions. Beside the enhancement of the branching ratios of 
class II processes, the inclusion of nonfactorizable contributions reduces the sensitivity to whether LL or NLL Wilson 
coefficients are used in the calculation. This is demonstrated in Fig. |3| by plottting the percentage difference between 
1^2 ^lii and \o,2^2\%ll as a function of eg,, taking e\ = 0. From the graph we see that a 10% contribution to eg 
reduces the difference between the branching ratios predicted by LL and NLL to about half that in naive factorization. 
Penguin diagrams contribute to processes of types b — > ccs and b — ► ccd. In both types, class II processes are 
affected only slightly by the penguin contributions. This is due to the destructive interference between the different 
terms, in the amplitude, generated by the penguins. On the other hand, for class I processes of type b — > ccs this 
cancellation does not happen and the decay amplitudes receive significant contribution from the penguin diagrams 



(see Table |HJ). 

As can be seen from Table [V] the best fit to the experimental data is produced by the Heavy-Quark Effective 
Theory (contained in the NRSX model), and lends support to the assumption of universality of the nonfactorizable 
contributions in B decays. 
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FIG. 1. The first few terms of 1/N C expansion, shown in a diagramatic language, that contribute to the decay amplitude of 

the process B — ■» D + tt~ . 
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FIG. 2. The annihilation diagram contributing to the decay amplitude of the process B — > D + tt 




FIG. 3. The Penguin Feynman diagrams for the process b — * ccs in the effective theory. 
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FIG. 4. (a) A contour plot of \ 2 m £i-£8 space using BSW II model, (b) A contour plot of \ 2 m £i-£s space using NRSX 
model. The lower graphs show a magnification of the region containing minimum 3 of x 2 i n the corresponding graphs. The 
inner closed curve represents a change of Ax = 1 from the minimum while the outer closed curve represents a change of 
A X 2 = 2. 
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FIG. 5. The percentage difference between |o2 £,2\ll an d l a 2 &\nll as a function of eg, taking ei — 
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TABLE I. The Wilson Coefficients in LL and NLL using the NDR. All values are evaluated at the scale jj, = 4.6 GeV, and 
for A|^ = 219 MeV. 



LL NLL 



d(p) 1.127 1.075 

C 2 (m) -0.286 -0.178 

C :i (p) 0.013 0.012 

C 4 (m) -0.029 -0.033 

C 5 (m) 0.008 0.009 

C 6 (m) -0.037 -0.039 

ai(/x) 1.032 1.016 

a 2 (/x) 0.090 0.180 

C 2 /ai -0.277 -0.175 

Ci/a 2 12.54 5.96 



TABLE II. Effective Wilson Coefficients in NLL, evaluated using A5_ = 219 MeV and the running quark masses at the 
6-quark mass scale. 

Cf 1.143 



Cf -0.322 



q — u q = c 

C? off 0.0184 + 0.0048 i 0.0197 + 0.0044 i 



'3 
1 
l 
J off 



C? eff -0.0407-0.0145i -0.0453 - 0.0132 i 



C| ett 0.0130 + 0.0045 i 0.0145 + 0.0044 i 

Cg cff -0.0522-0.0145^ -0.0568-0.0132^ 

ai 1.036 

a 2 0.059 

C 2 /ai -0.311 

Ci/a 2 19.37 
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TABLE III. AZ3f ac represents the percentage change in the branching ratio, assuming factorization, due to NLL values of 
the Wilson coefficients (column 4) and due to penguin diagrams (column 5). The last column represents the total change. 



Processes 



Class 



Model 



ASf ac : Change in branching ratio 
NLL effect Penguin effect Total change 



.etc. 



Type b — > cud 

~B° -> D+n~ ...etc. 

B -» D°W° . 

B~ -> D°%~ 
B~ -> D°p- 
B~ -* D°a~ 
B~ -* D* a %- 
B- -» D* a p- 
B- -> £>*% 



I 

II 

III 

III 

III 

III 

III 

III 



BSWII 
BSWII 
BSWII 
BSWII 
BSWII 
BSWII 



0.8% 
-56.7% 
-4.4% 
-2.1% 
-1.1% 
-6.0% 
-3.3% 
-1.8% 



0.8% 
-56.7% 
-4.4% 
-2.1% 
-1.1% 
-6.0% 
-3.3% 
-1.8% 



Type b — ► ccs 

B -> DD S 

B -> DD* 3 

B -» D*D S 

B -> D*D S * 

B -> KJ/i> . . . etc. 



I 
I 
I 
I 
II 



BSWII 



0.8% 
0.8% 
0.8% 
0.8% 
-56.7% 



-27% 
-36.0% 

6.4% 
15.1% 

0.2% 



-26.2% 

-35.2% 

7.2% 

15.9% 

-56.5% 



Type b — > ccd 

B -* DD' 
B -+ DD*~ 
B -* D*D~ 
B -> D*D*~ 

B — * nj/tp . . . etc. 



I 
I 
I 
I 
II 



BSWII 



0.8% 
0.8% 
0.8% 
0.8% 
-56.7% 



-23% 

-25.7% 

5.9% 

12.9% 

0.4% 



-22.2% 

-24.9% 

6.7% 

13.7% 

-56.3% 



TABLE IV. Predictions of the Fo and F\ form factors for the process B — » DD S in the five models considered in this work. 



Model 




Fo(m 2 Ds ) 



Fi{m\ 



BSWI 

BSWII 
NRSX |3 
AW 
ISGW 



0.756 
0.756 
0.653 
0.889 
0.832 



0.767 
0.849 
0.725 
0.937 
0.900 



TABLE V. The estimated nonfactorization parameters ei and eg in different models. The last column shows the value of \ 2 
per degree of freedom which indicate the goodness of fit. 



Model 



£i 



£s 



X 2 /d 



BSWI 

bsw nj§,|ij 

NRSX (32 

AW 
ISGW 



-0.037 ± 0.024 
-0.040 ± 0.024 
0.057 ±0.027 
-0.250 ±0.019 
-0.071 ± 0.024 



0.139 ±0.006 
0.137 ±0.006 
0.132 ±0.006 
0.171 ±0.007 
0.230 ± 0.008 



5.8 

0.6 
0.4 
5.7 
2.3 
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TABLE VI. The branching ratios predicted for a number of Cabibbo-favored B decays in BSW II model. The values in the 
second column were calculated by taking ei(po) = £s(mo) — whereas the values in the third column were calculated by taking 
£1(^0) = —0.040 ± 0.024 and eg((J,o) = 0.137 ± 0.006. The last column represents the available experimental measurements. 



Process 



Fac. 



Nonfac. 



Branching Ratio x 10 3 
3.4 ±0.2 

8.4 ±0.5 

9.5 ±1.0 

2.6 ±0.2 
7.6 ±0.4 
10.9 ±1.2 



Exp. B 



B" 
B° 
B° 
B° 
B° 
B° 



D+iT- 
D+p- 
D+a~ 
ZT+tT 
D* + p- 



4.1 ±0.1 
9.9 ±0.3 
11.3 ±1.0 

3.2 ±0.1 
9.0 ±0.2 
12.9 ±1.2 



Branching Ratio xlO 4 
0.77 ±0.16 
1.09 ±0.23 
0.43 ±0.09 
0.59 ±0.13 
0.13 ±0.03 
0.18 ±0.04 
0.53 ±0.11 
1.13 ±0.24 
0.52 ±0.11 
1.12 ±0.24 



3.0 ±0.4 
7.9 ±1.4 
6.0 ±3.3 

2.76 ±0.21 
6.7 ±3.3 

13.0 ±2.7 



B° 
B° 
B° 
B° 
B° 
B° 
B° 
B° 
B° 
B° 



L) 7T 

D*V 

D° v 

D*% 

D\> 

D*\' 

D°p° 

D*°p° 

D°u 

D*°w 



0.06 ±0.01 
0.08 ±0.02 
0.03 ±0.01 
0.05 ±0.01 
0.01 ±0.01 
0.01 ±0.01 
0.04 ±0.01 
0.09 ±0.02 
0.04 ±0.01 
0.09 ± 0.02 



< 1.2 
<4.4 

< 1.3 

< 2.6 
<9.4 

< 14 
<3.9 

< 5.6 
<5.1 

< 7.4 



B' 

B- 
B- 
B' 

B- 
B- 



D°p- 

D°a~ 

D*°n~ 

D*°p- 

D*°a7 



4.8 ±0.1 
11.1 ±0.3 
12.4 ±1.1 
3.8 ±0.1 

10.3 ±0.3 

14.4 ± 0.4 



Branching Ratio xlO 
5.4 ±0.3 
10.9 ±0.7 
11.6 ±1.2 
4.6 ±0.3 
10.9 ±0.3 
13.9 ±0.4 



Branching Ratio x 10 3 
8.5 ±1.9 
6.5 ±1.4 
7.1 ±1.6 
24.9 ± 5.2 



5.3 ±0.5 

13.4 ±1.8 

4.6 ±0.4 

15.5 ±3.1 



B -> DD 3 
B -> DD* S 
B -»• D*D S 
B -» D*D* 



10.3 ±2.2 
7.9 ±1.6 
8.4 ±1.8 

29.4 ± 5.9 



Branching Ratio xlO 4 
9.4 ±0.9 

5.0 ±0.6 
15.2 ±1.5 

9.1 ±1.1 



9.8 ±2.4 
9.4 ±3.1 
10.4 ±2.8 
22.3 ±5.7 



B -> KJ/ip 

B -+ Ki}j(2S) 
B -» K*J/iP 
B -» K*ii{2S) 



0.72 ± 0.05 
0.38 ±0.04 
1.17 ±0.09 
0.70 ± 0.07 



9.5 ±0.8 
5.8 ±1.2 
14.6 ±1.4 

9.6 ±2.5 



1G 



TABLE VII. The branching ratios predicted for a number of Cabibbo-favored B s decays in BSW II model. The values in the 
second column were calculated by taking ei(po) = £s(mo) — whereas the values in the third column were calculated by taking 
£i((M)) = —0.040 ± 0.024 and eg((J,o) = 0.137 ± 0.006. The last column represents the available experimental measurements. 



Process 



Fac. 



Nonfac. 



Branching Ratio x 10 a 

3.0 ±0.2 
7.3 ± 0.5 

8.3 ±0.9 
2.2 ±0.2 

6.4 ±0.5 

9.1 ±1.0 



Exp. 



B a 

B s 
B s 

B, 
B s 
Bs 



D+n~ 
D+p- 

Dt + TT- 

D*+p- 
D*+g- 



3.6 ±0.2 

8.7 ±0.4 
9.9 ±1.0 
2.6 ±0.1 
7.6 ±0.4 
10.8 ±1.1 



Branching Ratio x 10 4 
1.0 ±0.2 

1.4 ±0.3 
0.7 ±0.2 

1.5 ±0.3 



< 130 



B s -» D°K° 
Bs -> D*°K° 
B s -> D°K*° 
B 3 -> D*°K*° 



0.08 ± 0.02 
0.11 ±0.02 
0.06 ±0.01 
0.12 ±0.02 



B s -► riJ/i> 
B s - »?V(2S) 

Bs -► v'J/^ 
B s — t?>(2S) 

B s -> #(2S) 



0.19 ±0.02 
0.10 ±0.01 
0.22 ±0.02 
0.10 ±0.01 
0.79 ± 0.07 
0.47 ± 0.05 



Branching Ratio x 10 
2.5 ±0.3 
1.4 ±0.2 
2.9 ±0.3 
1.3 ±0.2 
10.3 ±1.1 
6.2 ±0.8 



< 38 



9.3 ±3.3 



TABLE VIII. The branching ratios predicted for a number of Cabibbo-suppressed B decays in BSW II model. The values 
in the second column were calculated by taking £i(/Uo) — £s(mo) = whereas the values in the third column were calculated by 
taking £i(/Cio) — — 0.040±0.024 and es(^o) — 0.137±0.006. The last column represents the available experimental measurements. 



Process 



Fac. 



Nonfac. 



Exp. 



Branching Ratio xlO 4 
3.0 ±0.6 
2.9 ±0.6 
2.4 ±0.5 
8.3 ±1.7 



B -> D+D~ 
B° -> D+D*- 
B° -» D*+D~ 



B 



D*+D" 



3.6 ±0.7 

3.5 ±0.7 

2.7 ±0.6 
9.9 ±2.0 



6.2 ± 3.6 H 



B~ 
B' 
B- 
B~ 
B' 
B~ 



-k J/rp 

p-_J/i> 

a 1 J/ip 

ir~ip(2S) 

p-il>(2S) 

a1%jj{2S) 



0.35 ± 0.03 
0.51 ±0.04 
0.23 ±0.02 
0.21 ±0.02 
0.33 ±0.03 
0.11 ±0.01 



Branching Ratio xlO 

4.6 ±0.5 

6.7 ±0.7 
3.0 ±0.3 

2.8 ±0.3 

4.4 ±0.5 

1.5 ±0.2 



5.0 ±1.5 
< 77 
< 120 
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